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■ Abstract. Big Bang nucleosynthesis (BBN) is the earliest sensitive probe of the 

values of many fundamental particle physics parameters. We have found the leading 
(— I I linear dependences of primordial abundances on all relevant parameters of the standard 

pi !■ BBN code, including binding energies and nuclear reaction rates. This enables us 

Q I to set limits on possible variations of fundamental parameters. We find that ''Li is 

i-^ ' expected to be significantly more sensitive than other species to many fundamental 

C/3 , parameters, a result which also holds for variations of coupling strengths in grand 

^ ' unified (GUT) models. Our work also indicates which areas of nuclear theory need 

further development if the values of "constants" are to be more accurately probed. 

The constancy over space and time of the coupling strengths and particle masses in the 
00 ■ Standard Model of particle physics is an assumption that should be tested [1]. Variations 
^ ■ may arise in a relativistically covariant theory due to the coupling of Standard Model 



particles to a scalar field whose cosmological value depends on time [2]. The possibility 
of time-varying couplings has been discussed since the first suggestions of dynamical 
dark energy [3j due to a scalar field evolving over recent cosmological epochs. 
■ Further motivation comes from possible signals of nonzero variation at redshifts 

' 0.5-4, arising from absorption spectra which probe the fine structure constant a and 
^ the proton-electron mass ratio = nip /me- The observational situation is contradictory, 

with both nonzero [H [5] and null [6] results, and debate about methods and statistical 
and systematic errors [7]. Bounds at different redshifts, and in different environments, 
probe possible variations without relying on particular models. We consider variations 
at the epoch of BBN {z ~ 10^'^), which is currently the earliest time at which theories 
of nuclear and particle physics can be compared to astrophysical observation. 

The study of particle physics via BBN faces two theoretical challenges. First, 
there are many variable parameters, to be compared with the small number of reliably 
observable abundances. Second, the uncertainty in how QCD parameters, in particular 
quark masses, affect nuclear forces, and thus nuclear binding energies and cross-sections. 
A derivation of nuclear forces from first principles, apart from the long-range attraction 
attributable to pion exchange, is lacking. Instead, various types of effective theories 
have been used which fit measured nuclear properties. 
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We present a systematic approach to these challenges. Firstly, each variation in 
particle physics or "fundamental" parameters is considered independently. Thus at 
linear order we can allow for any theoretical scenario in which such variations are subject 
to some unified relation. Second, we identify which nuclear properties and reactions 
have significant influence on the variation of the primordial abundances. This is done 
by varying every relevant nuclear binding energy and cross-section in a modified reaction 
integration code and noting the dependence of output abundances. This response to 
"nuclear" parameters is then to be related to the variation of fundamental parameters 
via nuclear theory. The connection between the different levels of physical understanding 
arises by simple matrix multiplication of "response matrices" . 

The dependence of primordial abundances with a = (D, '^He, ^He, ^Li, ^Li) on 
the variation of a set of nuclear physics parameters Xi is given by a response matrix C 
with matrix elements [8] 
_ dlnY, 
'''''' din X,- 

extracted by considering small variations of the Xi within the BBN code. This includes 
variation of the reaction rates which have a physical dependence on Xj. Then, the 
nuclear physics parameters X^ are related to a set of Standard Model parameters Gk 
via a second response matrix F with entries 
_ d\nX, 

This step requires, at present, theoretical assumptions, and contains substantial 
uncertainties in the dependence of nuclear binding energies on quark masses, which 
we discuss following Equations [7] and [81 The variation of abundances with respect to 
fundamental parameters Gk is then given by a matrix R, obtained simply via R = CF, 
with elements Vak'- 

AF, AGk 

= Tak-^. (3) 

We use mass-energy units where the QCD invariant scale Ac is kept constant. This 
is convenient for nuclear reactions where mass scales are determined mainly by strong 
interactions. Variations of dimensionful parameters are measured relative to Ac, for 
example considering the electron mass irtf. we implement a variation of nif./ Kf.- 
We vary with respect to the following thirteen "nuclear" parameters Xi. 

• Gravitational constant Gn 

• Neutron lifetime Tn 

• Fine structure constant a 

• Electron mass mg 

• Average nucleon mass m^v = {frtn + ni^jl 

• Neutron-proton mass difference Qat = '^n ~ 

• Binding energies of D, T, ^He, ^He, ®Li, ^Li, ^Be. 
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Table 1. Response matrix C, dependence of abundances on nuclear parameters. 



dlnYa/dlnXi 


D 


3He 


^He 


6Li 


^Li 


Gn 


0.94 


0.33 


0.36 


1.4 


-0.72 


a 


2.3 


0.79 


0.00 


4.6 


-8.1 




0.41 


0.15 


0.73 


1.4 


0.43 




-0.16 


-0.02 


-0.71 


-1.1 


-0.82 


Qn 


0.83 


0.31 


1.55 


2.9 


1.00 


niN 


3.5 


0.11 


-0.07 


2.0 


-12 


5d 


-2.8 


-2.1 


0.68 


-6.8 


8.8 


fix 


-0.22 


-1.4 





-0.20 


-2.5 




-2.1 


3.0 





-3.1 


-9.5 


-B4HC 


-0.01 


-0.57 





-59 


-57 













69 





Bju 














-6.9 


BjBc 














81 



Our results for these parameters are shown in Table [T], where the rows constitute the 
transposed nuclear response matrix C^; detailed discussion can be found in [S]. The 
largest sensitivity of abundances to nuclear parameters involves the variation of ^He, ^Li 
and ''Be binding energies. For instance the rate for ^He(a, 7)^Be with Q-value 1.59 MeV 
is very sensitive to changes in these (numerically large) binding energies. Q-values affect 
abundances both by influencing the reverse reaction rate as 
(0-^7)34^12 



oc e 



-Q/T 



(au) 12^34 

and via the kinematic dependence of cross-sections, 
reactions we have 



(4) 

For dipole radiative capture 



(5) 



aiE) oc E-^r^{Q + Ef 

whereas for 2 — 2 inelastic scattering or transfer reactions the dependence is 

a{E) oc ~ (g + Efl^ (6) 

where (3 is the outgoing channel velocity. For resonances, we scale their contributions 
to thermal reaction rates, which vary as e"^''/'^, by the appropriate power of [Q + Er). 

The dependence of reaction matrix elements on binding energies and on Q is in 
general not clear due to the lack of systematic effective theory. For the npd'') reaction 
we use the nuclear effective theory result of [lO] where dependence on Bj^ is explicit. In 
order to diagnose which reaction rates are important in the variation of final abundances, 
we varied each thermal averaged cross-section {av) by a temperature-independent factor, 
preserving the relation between forward and reverse rates. Apart from npd'^ and n p, 
only seven reaction rates strongly influence final abundances (discounting ^Li). In every 
case the dependences d\nYa/d\n{av)i are order unity or smaller [9]. Also, the ''He 
fraction does not depend on any reaction cross-section apart from n ^ p. 
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Table 2. Response matrix F, dependence of nuclear parameters Xi on fundamental 
parameters Gk- 



d\nXi/d\n Gk 




a 




me 




m 




1 

















a 





1 














Tn 





3.86 


4 


1.52 


-10.4 
















1 








Qn 





-0.59 








1.59 





mjv 

















0.048 


Bb 





-0.0081 











-4 


B'Y 





-0.0047 











-2.1/t 


BsHe 





-0.093 











— 2.3/3He 


BaRc 





-0.0304 











-0.94/4He 


Bqu 





-0.0541 











-1.4/6Li 


Bru 





-0.0459 











-1.4/7Li 


BjBc 





-0.0885 











-1.4/7Be 



Next we connect these nuclear parameters to fundamental parameters Gk at a higher 
energy scale. We consider the following six fundamental parameters: the gravitational 
constant Gn', the fine structure constant a; the electron mass me, the light quark mass 
difference 6g = m^ — m^, the averaged light quark mass m = (m^ + m„)/2 oc m^; and 
the Higgs v.e.v. (0). These are to be varied independently. The strange quark mass 
is omitted from our list because the present theoretical uncertainties in its influence on 
nuclear parameters are too high. 

Linear variations in the nuclear parameters Xi are then encoded in the matrix F 

defined in ([2]): our estimates of F [9] are shown in Table [2l We estimate the leading 

dependence of nuclear binding energies on the pion mass, due to the pion-mediated 

contribution to the nucleon-nucleon potential, by scaling up the dependence of the 

deuteron binding energy. For the deuteron, 

r 

A In Bd = r A In = -A In m (7) 
with r ^ —8 ± 2 [llj. Then for a nucleus with mass number Ai we estimate 



1^ = MAi - l)^r ^ -O.miAi - 1) 
omTT m-n 



(8) 



where the constants are expected to be of order unity, but may vary in magnitude by 
factors of a few. This is our main source of theoretical uncertainty, in the case where m 
varies significantly. However, we expect our main finding of a much larger dependence of 
^Li than of D and ^He to be robust, barring fine-tuned cancellations. Table [3]then shows 
the dependences of abundances on fundamental parameters, encoded in the matrix R. 
For the m-dependences we have given the values which arise when setting all to unity. 

We may now set bounds on the variation of each fundamental parameter, considered 
in isolation. Out of three observational determinations of primordial abundances (D, 
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Table 3. Response matrix R, dependence of abundances Yi on fundamental 
parameters Gk- 



dinYa/dinGk 


U 


3tt 


He 


6t • 

°Li 


7t • 

'Li 


Gn 


0.94 


0.33 


0.36 


1.4 


-0.72 


a 


3.6 


0.95 


1.9 


6.6 


-11 


(0) 


1.6 


0.60 


2.9 


5.5 


1.7 


me 


0.46 


0.21 


0.40 


0.97 


-0.17 




-2.9 


-1.1 


-5.1 


-9.7 


-2.9 


171 


17 


5.0 


-2.7 


-6 


-61 


7] 


-1.6 


-0.57 


0.04 


-1.5 


2.1 



Table 4. Allowed individual variations of fundamental couplings. 



-19% 


< 


A In Gn 


< 


+ 10% 


-3.6% 


< 


A In a 


< 


+1.9% 


-2.3% 


< 


Aln(0) 


< 


+1.2% 


-17% 


< 


A In nie 


< 


+9.0% 


-0.7% 


< 


A\n6q 


< 


+1.3% 


-1.3% 


< 


A ln?fi 


< 


+1.7% 



^He and ^Li) the observed ''Li abundance is a factor two to three smaller than standard 
BBN theory predicts. Thus we use D and ^He for the purpose of constraining allowed 
variations. For deuterium we take 2a limits from [12j; for ^He we consider instead the 
"conservative allowable range" of [13]. The resulting constraints are given in Table HI 

We also consider unified scenarios where the variations of fundamental couplings 
satisfy relations that reduce the number of free parameters. If there is a single underlying 
degree of freedom which varies, written as a dimensionless scalar cp, then variations of 
Gk can then be written as a vector: AlnG^ = dkAip, where dk are a set of numbers 
characterising a particular unified model. We then obtain 

Alny; = (CF),fc4A<^. (9) 

One may eliminate A(f in favour of the variation of another parameter, such as a. 
We considered a grand unified theory with unified gauge coupling ax, broken at 
the scale Mx to the Standard Model group. The observable couplings of QCD and 
electromagnetism are related to ax via renormalization group flow and electroweak sym- 
metry-breaking. We take for simplicity the Planck mass fixed relative to the unification 
scale, A(Mp/Mx) = 0, and the Yukawa couplings to be constant, thus electron and 
quark masses are proportional to (0). The variation of (0) is parameterised as 



(0) 



const- 



/ Ac 



(10) 



Mx " \Mx, 

We found two cases of interest. In the first, the Higgs v.e.v. is proportional to the 
unification scale, thus 7 = 0. Then we find [9j 

Aln(GN,a, (0),me,5g,m) ~ (64.5, 1, -32, -32, -32, -32)A Ina. (11) 
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We then obtain variations of abundances 

A\n{Yj,,Y3ne,Yp,Yeu,Y7u) ^ (-450, -130, 170, 380, 1960) A Ina. (12) 
Note the strong variation of ''Li compared to other abundances. Here, a fractional 
variation of a by about —2.5 x 10^^ would bring abundances within 2a observational 
bounds. Secondly, we consider a case when the Higgs v.e.v. and fermion masses vary 
more rapidly (with respect to Mx) than the QCD scale: thus 7 > 1. We take 7 = 1.5 
and find variations of fundamental couplings 

Aln(GN,a, (0),me,5g,m) ~ (87, 1, 21.5, 21.5, 21.5, 21. 5)A Ina. (13) 

The variations of abundances are then 

A\n{Yj,,Ysne,yp,Yeu,Y,u) ^ (430, 130, -65, -60, -1420)A Ina. (14) 

In this scenario a fractional variation of a by about +4 x 10"'^ brings theory and 
observation into agreement within Icr bounds. Thus the "lithium problem" may be 
considerably ameliorated by allowing variable "constants" . 

In both scenarios, as ^Li decreases, "^He also decreases, but the D abundance 
increases. Then we would expect the true primordial ^He fraction to be at the lower 
side of its currently allowed range, while D is near its upper observational limit. 

In summary, primordial nucleosynthesis is a unique window on the very early 
universe, setting the earliest or most stringent bounds on many effects in cosmology 
and particle physics. The possibility that particle masses and couplings could have 
a space-time dependence can thus be tested, and we present a systematic treatment 
of potential signals and bounds at BBN. In general there are more variable parameters 
than observables, but if unification of gauge couplings is assumed, definite limits emerge, 
and preferred models may be found where the predicted ''Li abundance is strongly 
suppressed, while the helium fraction decreases and the deuterium abundance increases, 
within current observational limits. 
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